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Abstract. A brief review of the phenomenological studies in the field of heavy ion collisions 
based on the Color Glass Condensate theory and, in particular, of those relying in the use of 
the BK equation including running coupling effects is presented. 



1. Introduction 

The main goal of the experimental programs on ultra-relativistic heavy ion collisions at RHIC 
and the LHC is the production and characterization of the Quark Gluon Plasma (QGP). A 
crucial ingredient for the proper understanding of such collisions is the knowledge of the structure 
of the wave functions of the colliding nuclei, typically encoded in their parton distributions. At 
high energies, nuclei do not behave as a mere incoherent superposition of their constituents 
nucleons. Rather, coherence effects are important and modify not only their partonic content, 
but also the underlying dynamics of particle production in scattering processes. A careful 
distinction of such initial state effects from those originating from the presence of a QGP, or final 
state effects, is of vital importance for a proper characterization of the matter produced in heavy 
ion collisions, as they may lead sometimes to qualitatively similar phenomena in observables of 
interest. 

The Color Glass Condensate effective theory provides a consistent framework to study QCD 
scattering at high energies (for a review see e.g. [HE])- The main physical ingredient in 
the CGC is the inclusion of unitarity effects through the proper consideration of non-linear 
recombination effects, both at the level of particle production and also in the quantum evolution 
of hadronic wave functions. Such effects are expected to be relevant when nuclei (or hadrons, in 
full generality) are proven at small enough values of Bjorken-x. In that regime gluon occupation 
numbers are very large and gluon self-interactions become highly probable, thus taming, or 
saturating, further growth of the gluon densities. While the need for unitarity effects comprised 
in the CGC is, at a theoretical level, clear, the real challenge from a phenomenological point of 
view is to assess to what extent they are present in available data. Such is a difficult task, since 
different physical mechanisms concur in data, and also because the limit of asymptotically high 
energy in which the CGC formalism is developed may not be realized in current experiments. 
In that sense, the calculation of higher order corrections to the CGC formalism has supposed 
important leap forward in sharpening the CGC as an useful phenomenological tool. 

2. Some recent developments in the CGC: Running coupling BK equation. 

The leading order BK-JIMWLK equations resums soft gluon emission in the leading logarithmic 
(LL) approximation in a s \nl/x to all orders, besides of including non-linear terms required by 



unitarity. At such degree of accuracy, the theory is incompatible with data. Indeed, the in- 
dependence of the saturation scale emerging from the LL equations is Q 2 s {x) « Qq(xq / x) xLL , 
with \ LL ~ 4.8a s [3], while phenomenological studies of electron-proton scattering and heavy ion 
collisions indicate that an evolution speed A ~ 0.2 0.3 is preferred by data. Such insufficiency 
of the theory has been partially fixed by the calculation of running coupling corrections to the 
BK-JIMWLK equations through the inclusion of quark loops to all orders [HE]. Among other 
interesting dynamical effects, running coupling effects tame the growth of the saturation scale 
down to values compatible with experimental data [6]. Due to the complexity of the JIMWLK 
equations, in phenomenological works it is more feasible to solve the BK equation, more tractable 
numerically, which corresponds to their large- N c limit. It reads 

dM Q Y Y) = J d 2 ri^ mn (r,r 1 ,r 2 )[AA(r 1 ,y)+AA(r 2 ,y)-AA(r,y)-AA(r 1 ,y)AA(r 2 ,y)] , (1) 

where M(r, Y) is the dipole scattering amplitude on a dense target, Y = lnxo/x the rapidity 
and r the dipole transverse size. It turns out that running coupling effects can be incorporated 
to the evolution equation through just a modification of the evolution kernel, referred to as 
K run in Eq. (TjQ). As discussed in [6], there are several ways of doing so, since the separation of 
running coupling and conformal terms is not uniquely defined. The evolution kernel proposed in 
[5] minimizes the role of higher order corrections, making it better suited for phenomenological 
applications (see [6] for an extended discussion on the subject). Finally, Eq. (P) needs to be 
suplemented with initial conditions, which can be choosen to be of the McLerran-Venugopalan 
type [7]. This introduces two free parameters: The value xq where the evolution starts and the 
initial saturation scale Qq. Finally, the unintegrated gluon distribution entering the different 
production processes discussed below is related to the dipole amplitude in Eq. (TjQ) through a 
Fourier transform (see Eq. ([5])). Besides of the works in heavy ions described below, the rcBK 
equation has been shown to successfully account for the x-dependence of structure functions 
measured in electron-proton scattering [8J. 

3. Total multiplicities in heavy ion collusions 

The smaller than expected total multiplicities observed at RHIC suggest strong coherence effects 
in the mechanism for particle production in heavy ion collisions at high energies. The CGC offers 
a natural explanation of this observation: The total flux of scattering centers (gluons) entering 
the collisions is significantly reduced due to saturation effects in the wave function of the colliding 
nuclei. Such idea is realized in the phenomenological KLN model [S], which relies in the use 
of fct-factorization. There the total multiplicities in central collisions rise proportional to the 
saturation scale of the colliding nuclei, dN/drj oc Q 2 sA - Even though the use of ^-factorization 
in A+A collisions is not justified, the good description of the energy, rapidity and centrality 
of multiplicity data yielded by the model lends support to the underlying physical picture. An 
important improvement of the KLN model was presented in [TO], where the phenomenological 
model for the unintegrated gluon distributions was replaced by exact solutions of the running 
ccoupling BK equation, also getting a very good description of RHIC data, as seen in Fig. 1. 
Besides of reducing the degree of modeling needed to describe data, the use of well controlled 
QCD tools as the rcBK equation allows for a better controlled extrapolation to LHC energies. 
Taking the centrality class — 6% at RHIC as a baseline, the following mid-rapidity multiplicities 
at the two expected LHC energies in most central Pb+Pb collisions were obtained in [10J: 
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The error bands in Eq. ([2]) originate mostly from the uncertainties in the initial conditions for 
the rcBK equation. 




Figure 1. Charged particle multiplicity densities in Au+Au collisions at RHIC (data from 
[TT] ) and predictions for Pb+Pb collisions at ^/s^N = 2.76TeV (left). Mid-rapidity charged 
multiplicity as a function of the collision energy (left). 



4. Single inclusive particle production and nuclear modification factors 

Particle spectra in A+A collisions are expected to be strongly modified due to the presence of 
a QGP. Thus, proton (deuteron) nucleus collisions, where no QGP formation is expected, are 
better suited for the exploration of initial state effects. For any production process, nuclear 
effects are typically evaluated in terms of the nuclear modification factors: 
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where N co u is the number of collisions. If high-energy nuclear reactions were a mere incoherent 
superposition of nucleon-nucleon collisions, then the observed R p a should be equal to unity. 
However, RHIC measurements in d+Au collisions (or peripheral Au+Au collisions) |12t [T3] 
feature two distinct regimes: At mid rapidities the nuclear modification factors exhibit an 
enhancement in particle production at intermediate momenta pt ~ 2-^4 GeV. In turn, a 
suppression at smaller momenta is observed. However, at more forward rapidities such Cronin 
enhancement disappears, turning into an almost homogeneous suppression for all the measured 
values of pt- According to 2 — > 1 kinematics, the x-values probed in the projectile and target are 
x i(2) = ( m t /yfs) ■ Thus, for RHIC energies, mid rapidity implies moderate x ~ 0.1 -j- 0.01, 
probably not small enough to justify the applicability of the CGC. Indeed, mid-rapidity data has 
been analyzed through different formalisms and techniques, from collinear perturbation theory 
|14| to independent multiple scatterings [15] or CGC approaches [16] . all of them reaching a 
reasonably good agreement with data. Data collected in the deuteron fragmentation region offer 
a cleaner opportunity to explore CGC effects, since the nucleus is proben at muchsmaller values 
of x. There, the CGC formulation of single particle production takes on a relatively simple form 
[IT]: 

dN h K \- f l dz r 2 ~ / p t \ 2 



+ xifg/ p (x 1 ,pf)N A (»2)— J D h/g (z,pt) , (4) 



where pt and are the transverse momentum and rapidity of the produced hadron, and fj_/ p 
and Dfon refer to the parton distribution function of the incoming proton and to the final-state 
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Figure 2. Negative charged hadron and neutral pions at forward rapidities measured by the 
BRAHMS [12] and STAR |13j Collaborations in p+p (left) and d+Au minimum bias (right) 
collisions compared to our calculation [15] . 

hadron fragmentation function respectively. Thus, in the forward region the projectile is in the 
dilute regime and characterized by its parton distribution functions, while the nucleus is deep 
in the saturation region and characterized by unintegrated gluon distributions taken from the 
solutions of the rcBK equation: 



where k refers to transverse momentum. With this set up we reach a very good description 
of forward neutral pions and negatively charged hadrons yields as measured by the STAR and 
BRAHMS Collaborations respectively in d+Au minimum bias and in p+p collisions, as shown 
in Fig 2. All the details of the calculation and fit parameters can be found in [18 1. 

By simply taking the ratios of the corresponding spectra, we get a very good description of 
the nuclear modification factors at forward rapidities. It should be noted that we use the same 
normalization as the experimentalist do in their analyses of minimum bias d+Au collisions, i.e. 
we fix N co u = 7.2. Physically, the observed suppression is due to the relative enhancement of 
non-linear terms in the small- x evolution of the nuclear wave function with respect to that of 
a proton. However, it has been argued that the observed suppression at forward rapidities is 
not an effect associated to the small values proven in the nuclear wave function but, rather, 
to energy- momentum conservation corrections relevant for xp — > 1 Such corrections are not 
present in the CGC, built upon the eikonal approximation. Thus, the energy degradation of 
the projectile parton through either elastic scattering or induced gluon brehmstralung would be 
larger in a nucleus than in proton on account of the stronger color fields of the former, resulting 
in the relative suppression observed in data. A successful description of forward ratios based on 
the energy loss calculation was presented in [19J. 

Thus, at the level of only this observable, one is in an inconclusive situation. Disentangling the 
underlying dynamics of the suppression requires either data at larger energies or, alternatively, 
the study of more exclusive observables. The huge leap forward in collision energy reached at the 
LHC allows for an exploration of small- x effects already at mid-rapidity. There, both the target 
and projectile x are small, and energy loss effects associated to large-x^ effects are expected to be 
small. In Fig 4. we present our CGC predictions for the nuclear modification factor for negative 
charged hadrons in p+Pb collisions at two LHC energies. Our curves correspond to rapidities 
y = 2 and larger. Technical difficulties related to the intrinsic asymmetry of the formalism used 
for particle production prevent us from calculating the ratios at mid-rapidity. However, the 
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Figure 3. Nuclear modification factors at forward rapidities in minimum bias d+Au collisions 
in the CGC [18]. 



0.7 r 

0.6 r 



;.8 TeV y h =2, 4, 6 
N onll =3.6 



. hybrid formalism 
kt- factorization 



P T (GeV) 



\(s NN =6.2TeV y h =2,4,6 
N mll =3.6 



. hybrid formalism 
- kt-factorizatiori 



P T (GeV) 



Figure 4. CGC predictions for the nuclear modification factor in p+Pb collisions at two LHC 
energies and rapidities yh = 2, 4 and 6. 



smooth rapidity dependence suggests that a large suppression ~ 0.6 is also expected at mid- 
rapidity in the LHC. It should be taken into account that the normalization taken to produce 
the curves in Fig was N co u = 3.6. 

5. Forward di-hadron correlations 

As said before, disentangling CGC from other dynamical effects demands the analysis of more 
exclusive observables. The opportunity is provided by the recent measurement of forward di- 
pions correlations in d+Au collisions by the STAR collaboration [20] . The experimental data 
for the coincidence probability CP(Acf)) feature a clear suppression of the away-side correlations 
with respect to the ones observed p+p collisions. This suppression is naturally explained in the 
CGC: The two pions experimentally observed originate from a valence quark-gluon system in 
perfect back-to-back correlation. The quark-gluon system is put on-shell through the interaction 
with the nucleus, as a result of which the quark and gluon also acquire a transverse momentum 
of the order of the saturation scale of the nucleus. When that scale, which marks the onset 
of non-linear effects, is comparable to the initial transverse momenta of the quark and gluon, 
their intrinsic angular correlation is washed out. Finally, the outgoing quark and gluon fragment 
independently into hadrons. A recent calculation CGC [21] provides a very good comparison 
with available data. There, the required information on the nuclear wave function is built, 
through the use of the gaussian approximation, from the two-point function Eq. ([5]) constrained 
by the analyses of single inclusive spectra, making it a parameter-free calculation. This study 
provides, arguably, the most compelling evidence for the presence of CGC effects in available 
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Figure 5. Comparison of the CGC calculation of |21j with data for the coincidence probability 
at forward rapidities in d+Au (blue) and p+p collisions (red). 
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